Fock space expansion of cr meson in leading-A^c 
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Abstract 

We examine the leading-A^c behavior of the masses and transition matrix elements of some low-lying, few-particle 
configurations in QCD. A truncation of the Fock space produces an effective, symmetric Hamiltonian that we diago- 
nalize. The lowest eigenvalue is identified as the cr meson if the Hamiltonian is chosen to represent the scalar sector 
As an application, the coefficients of the Nc powers are then fit to two-loop Unitarized SU(2) Chiral Perturbation 
Theory results for the cr mass and width as a function of the number of colors, and we show that those results can 
be accommodated using the QCD A^^. dependence previously derived for matrix elements, without the need for un- 
natural parameters or fine tunings. Finally, we show a very preliminary good quality fit, estimating the proportion of 
tetraquark/molecule-like (dominant), ^^-like (subdominant) and exotic-like (marginal) configurations in the cr. 
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1. Motivation 

A low-mass cr meson was introduced in 1955 111] as 
an auxiliary device that has turned out to be very use- 
ful to explain the intermediate-distance attractive part in 
the nuclear potential. For long time debated, the mass 
and width of this meson (decaying almost always to nn, 
the only open strong decay channel) have been recently 
pinned down with very good precision employing differ- 
ent methods. The results of these analysis are in agree- 
ment with each other and some are shown in table [1] 

Table 1 : Precise determinations of the cr meson mass and half-width (in MeV). 



M 


r/2 


Refs. 


452(12) 


260(15) 


[2] 


458(15) 


262(15) 


[3] 


441(12) 


272(12) 


[4] 



It behooves one to understand the composition of this 
meson in terms of the fundamental QCD degrees of 
freedom, quarks and gluons. Of current interest is the 
decomposition of states in terms of a Fock space expan- 
sion Js] and this we address in the present brief report. 
This expansion reads 

|o-> = ^ J" {aqqlqq) + Ogglgg) + aqqqq\qqqq) . . .) (1) 



where the sum/integral signs remind us of spin, momen- 
tum and other degrees of freedom that, for simplicity, 
we will further omit in our notation; that is, we consider 
the summed amplitude over each Fock subspace 



k) = Ogijlqq) + agglgg) + agggglqqqq) . 



(2) 



'Speaker 



This expansion in terms of quarks and transverse glu- 
ons is well defined in Coulomb gauge QCD that 
can be formulated without ghosts nor longitudinal glu- 
ons. At least for heavy mesons decaying to open- 
flavor channels, the intrinsic qq component can be iden- 
tified in a model-independent way ff\. The setback 
of this full quantum-mechanical answer is that it is 
frame-dependent, presumably defined in the rest frame 
of the hadron fs']. This makes it less attractive for light 
hadrons where speeds can be large. 

The 1 /Nc expansion around Nc = 3 offers more Um- 
ited information: it can only separate classes of equiv- 
alence of states whose mass and decays behave in the 
same way under A^^, but the information obtained is 
useful also for light quarks. Thus our states in Eq.(|2|i 
should be understood as qq-liks, gg-like, etc... although 
for simplicity we are calling them qq, gg, etc... 

2. Matrix elements in leading-A^c 

We now consider what configurations may play 
an important role for the cr meson (and other light 
scalar mesons). First, the cr is a very broad reso- 
nance in nn scattering, deep in the complex plane. 
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so one can describe part of its nature as a pion-pion 
correlation (distortion of the density of states, or for 
brevity, 'molecule'), or equivalently as the leading- 
Nc color analysis is concerned, a tetraquark. Thus, 
although we denote it by qqqq, we actually mean 
"tetraquark/molecule". 

Next, at some level one expects to find (and indeed 
finds as will be shown in figure |2] below) a qq compo- 
nent, that would correspond to the 1 Ge V quark model's 
^Pq configuration. 

Finally, one might expect more exotic configurations 
such as glueballs or baryonium-type multiquark corre- 
lations to also play a role. 

We give the leading A^^. behavior, N^, of the mass and 
7T7T width of these various configuration in table |2] 

Table 2: Leading-W^ scaling of mass and width of various QCD Fock-states. 



State 


M 


r 

^ 7T7T 


nn, qqqq 


0(1) 


0(1) 




0(1) 


0(1 /A',) 


gg 


0(1) 


0(1 /N^) 


{N, - l){qq) 


0(N,) 





We illustrate the color computation of one of these ma- 
trix elements (the glueball to two pion transition G 
tttt) in figure [1] 





Figure 1 : Left: the impulse diagram for the transition of a glueball to two 
pions already yields the leading-W^ behavior of the entire amplitude as shown 
by t'Hooft. Right: color flow of the same diagram. The line crossing reveals the 
1 /N^ suppression. 



A way to establish the counting (left diagram in the 
figure) is to observe that the color-singlet two-gluon 
wavefunction, properly normalized, is —f==. Each 



of the two vertices carry (this scaling of the color 
charge guarantees that higher-order diagrams scale in 
the same way under A^^)- Finally, the pion wavefunc- 
tions in the final state combine a quark and antiquark to 



the matrix element 



form a color singlet . 

The net result for 
tr{T"T")l{Nl ^JNf^), suppressed as IjNc. 
the width is proportional to the matrix element G — > 
squared, the corresponding entry in table|2]follows. 



IS 
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3. Effective Hamiltonian 

We take one state from each class qq, gg, qqqq^ to 
build a discrete 3x3 effective Hamiltonian H. We will 
not attempt to calculate this Hamiltonian from theory, 
but we factor the leading-A^c behavior of its matrix el- 
ements, known from a Nc analysis, and leave the pre- 
coefficients as free parameters. Diagonalization of H 
yields three eigenvalues. We identify the lowest one as 
the cr. 

To describe the width, our three model states have 
to be coupled to the pion-pion continuum. For this we 
employ a Feshbach decomposition 1, 17.1 in terms of a 
P subspace (our three discrete states) and a Q subspace 
(two free pions in an arbitrary relative momentum state). 
The full Hamiltonian in the total space 



H = 



Hpp 

Hr 



Hn 



PQ 



(3) 

iQP nQQ ■ 

is then restricted to the discrete P subspace via the resol- 
vent in 2-space with appropriate boundary conditions 

1 

HpQ— — — —Hqp 



H' 



■ff 



pp 



pp 



E-Hr 



(4) 



iQQ + le 

We actually do not need to calculate the integral over 
pion configurations in the rightmost term; all we need is 
to extract its leading-A^c behavior. The effective Hamil- 
tonian is finally a symmetric (because of CP invari- 
ance), non-Hermitian (because of the Fock-space re- 
striction) 3x3 complex matrix, that has therefore 12 free 
parameters in leading Nc (the exponents being known). 



(5) 



The diagonal yS are given in Table 2, but for brevity, all 
others - already calculated with the procedure explained 
above- will be given somewhere else. We express the 
lowest eigenvalue (now complex) in terms of these un- 
specified parameters hij. 

4, Inverse Amplitude Method 

We match our Nc results to a dispersive analysis, 
the well-established Inverse Amplitude Method (lAM) 
The method proceeds by writing a dispersion 
relation for G = '^-j^ where t is the scalar, isoscalar 
pion-pion scattering amplitude and t - fi'' + f^"*' -i- f*^^ . . . 
is its ChPT expansion. This reads 

G(s) = G(0) + G'{Q)s + ^G"(0)s^ 



r , ImG(i') 
— ds'— — — + LC(G) + PC(G) 
^ Jrc s'\s' - s) 



' For this analysis we discard the baryonium-like configuration, 
that will be reexamined in an upcoming work. 
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Figure 2: Mass (top) and width (bottom) of the tr meson in Unitarized Chiral 
Perturbation Theory (dots) as a function of the number of colors. The solid line 
is the fit to these data points in the three-state model that includes a inolecule- 
like state, and two intrinsic qq and gg states. Squares (puiple colored online) at 
A^^. - 3 represent the precise dispersive determinations from tablefTl 



in terms of a left Mandelstam cut (LC), pole contribu- 
tions (PC), due to the Adler zeroes of t and analyzed 
elsewhere ll l[| . polynomial subtractions to assist con- 
vergence, and a dispersive integral over the right, unitar- 
ity cut. The two first terms in the first line are approxi- 
mated in Chiral Perturbation Theory. Their contribution 
to the amplitude in the physical region with E > 2m„ is 
very small since they lay well to the left of this thresh- 
old in the complex .s-plane. The polynomic subtraction 
is represented exactly in the chiral expansion. 

This leaves the dispersive integral on the second line. 
The nice feature about it is that, in the elastic region 
for nn scattering, with 2mK > E > 2m„, the imaginary 
part is exactly known, Im G = -Im t4, and this formula 
provides a good approximation further up to £ ^ 1.2 
GeV. 

The dispersion relation can then be turned into simple 
algebraic expressions to the order desired, and while the 
approximation to order p"* has been much exploited, we 
here employ the order in the expansion, that yields 
(not writing down the pole contributions for simplicity) 



(f2 - f4 + tl/t2 - fe) 



(6) 



l[ix 10^) 


-5.4 


ri(x 10^) 


-0.6 


/^(x 10^) 


1.8 


riix 10^) 


1.5 


/^(x 10^) 


1.5 


r3(x 10^) 


-1.4 


F.ix 10^) 


9.0 


r4(x 10^) 


1.4 






rsix lO'*) 


2.4 






r6(x 10^) 


-0.6 



Table 3: Two-loop lAM LECs employed, corresponding to the fit p 
as qq in fT3ll . We rely on 5 1/(2) chiral perturbation theory, while it 
is known that the subthreshold coupling of the cr to the closed KK 
channel is large 01 . A pure color analysis is blind to these flavor 
details, and for the time being the KK component must be understood 
as included in the (dominant) meson-meson component analyzed. 



ing are thus simply obtained as zeroes of the denomi- 
nator, and the cr mass and width extracted thereof. The 
dependence Nc dependence of the resonance parame- 
ters then obtained |12]|13| by changing the chiral low 
energy parameters following their model independent 
ChPT/QCD description 114] (that follows from count- 
ing flavor traces). For this cr meson application we 



fn J 3' , li 



/,4- for / = 1 ...4 and 



need only f„ 

fi — > fi (t)" ^'-'^ ' ~ 1 ... 6. The values of the low en- 
ergy constants at A^^. = 3 are given in table[3] 

Figure |2] does not exceed A^^. ^ 20. The lAM is 
reliable near A^^. - 3 where the resummation of s- 
channel rescattering effected by the JAM is the domi- 
nant physics, and the unitarity cut dominates the dis- 
persion relation. Also a flavor singlet Goldstone boson 
is not necessary in the effective Lagrangian for modest 
A^e, as the axial anomaly is Nc suppressed. Our results 
should be understood as a leadin^-Nc expansion around 
Nc - 3, never as Nc — » °°, see OlSll . 



The poles (elastic resonances) in pion-pion scatter- 



5. Results 

The JAM results ifTill for varying-A^c are displayed 
in figure |2] The mass Mo- initially increases towards 
1 GeV with growing A^c, but then saturates to a more or 
less constant value after A^^ ^ 8. The half-width ro-/2 
also starts increasing with Nc, but around A^c ^ 6, it de- 
creases as it is expected for a qq meson. This behaviour 
has been interpreted ifisll as a signal of the onset of the 
intrinsic, qq subdominant component of the cr meson 
and seems to be needed to ensure fulfillment of local 
duality [16J . The behavior of the width at lower A^c is 
characteristic of a molecule or tetraquark component. 

It is plain from a comparison with table [T] that none 
of the intrinsic QCD states by itself can reproduce this 
behavior. One therefore needs to consider the mixing 
between different configurations. 

Employing our 3x3 effective Hamiltonian, and 
varying Nc in the known manner, we fit the free pre- 



3 



coefficients hij. We restrict the fit to the parameter sub- 
space yielding one light scalar only (in agreement with 
Unitarized Chiral Perturbation Theory), the other two 
being above 1 .2 GeV, and we make no further statement 
about them since at that point they escape the reach of 
our method. We further examine naturality in the A^^ 
pre-coefficients, in the sense that a coefficient of the 
known Nc powers is said to be natural if h/j € ( 1 /Nc, Nc). 

This equation guarantees that the naive Nc counting 
works for the effective Hamiltonian. It could acciden- 
tally happen that one of the pre-coefficients of a sub- 
leading term was very large and the lowest orders in the 
1 /Nc expansion were therefore not a reliable approxima- 
tion. Our naturality assumption means that we discard 
this case and believe the Nc counting as is. 

Our very preliminary best fit is shown in figure |3] 
where we plot the probabilities \a\^ of finding one of 
our three model states in the lowest eigenvector. 



1 
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Nc 



Figure 3; Fock .space decompo,sitioii of the cr meson in a discrete model 3- 
state subspace, imposing naturality for the effective 3x3 Hamiltonian, whose 
pre-coefiicients are obtained fitting leading-W^. to Unitarized Chiral Perturbation 
Theory. The dominant component behaves under N^- as a qqqq through A'^. =^ 6, 
then the subdominant ^^-like takes over for larger Nc . The glueball-like com- 
ponent stays always at or below the 10% level. 

The graph encodes all information that can be ex- 
tracted from the Nc counting and naturality alone with- 
out increasing the model space. Given the large number 
of parameters in the minimization, we would be cau- 
tious about extending P beyond, say, four states. 

It is rewarding that a good fit can be found (the soUd 
line in figure |2]i with such a simple model. 

One could also inquire whether alternative fits that 
make the exotic (glueball) component dominant in the 
cr expansion could be found as some authors indicate 
that this might be the case [2i- This is ongoing work 
that will be reported in a follow-up publication. 

A similar analysis could also be carried out for 
baryons and in fact, another problematic state in the low 
hadron spectrum, the A(1405), is being examined ifigll . 
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